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Graphene [If is a sturdy and chemically inert material exhibiting an ex- 
posed two-dimensional electron gas of high mobility. These combined proper- 
ties enable the design of graphene composites, based either on covalent [2] 
or non-covalent [3] coupling of adsorbates, or on multilayered structures [4]. 
These systems have shown novel tunable electronic properties such as bandgap- 
engineering [3], reversible metal/insulating transition [2], |4] or supramolecular 
spintronics [5j. Tunable superconductivity is expected as well [6J, but experimen- 
tal realization is lacking. Here, we show experiments based on metal/graphene 
hybrid composites, enabling the tunable proximity coupling of an array of super- 
conducting nanoparticles of tin onto a macroscopic graphene sheet. This material 
allows a full electrical control of the superconductivity down to a strongly insu- 
lating state at low temperature. The observed gate control of superconductivity 
results from the combination of a proximity-induced superconductivity genera- 
ted by the metallic nanoparticles array with the two-dimensional and tunable 
metallicity of graphene. The resulting hybrid material behaves, as a whole, like 
a granular superconductor showing universal transition threshold and localiza- 
tion of Cooper pairs in the insulating phase. This experiment sheds light on the 
emergence of superconductivity in inhomogeneous superconductors, and more 
generally, it demonstrates the potential of graphene as a versatile building-block 
for the realization of novel superconducting materials. 

Even though intrinsic superconductivity in doped [6] graphene has been proposed, it has 
not yet been experimentally shown. However graphene can carry supercurrent by means 
of proximity effect [7] , in which charge carriers in a non-superconducting material acquire 
superconducting correlations in the vicinity of a superconductor. In the present study, the 
proximity effect in graphene is not locally generated from the contacting electrodes, like it 
is done in most mesoscale experiments [7] , but rather by coupling its surface to a macrosco- 
pic 2D network of superconducting clusters [8|. As opposed to covalent functionalization of 
graphene, which strongly affects the density of state, non-covalent coupling of adsorbates on 
graphene is useful for designing materials with novel electronic functions, since the graphene 
keeps most of its exceptional electronic properties, while gaining others coming from the 
coupling elements (Fig.l). Previous experiments involving decoration of metal clusters on 
exfoliated graphene have shown [9] that such a hybrid system can exhibit a gate-tunable 
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Berezinsky-Kosterlitz-Thouless (BKT) transition towards a fully two-dimensional supercon- 
ducting state with a critical temperature related to the normal state resistivity. In that expe- 
riment, the mean free path and the superconducting coherence length exceeded the average 
gap width separating neighboring tin nanoparticles, resulting in homogeneous 2D super- 
conductivity. Here, the main and crucial difference is the use of centimeter-scale graphene 
layers grown by chemical vapor deposition (CVD) instead of previously used micron-scaled 
exfoliated graphene. 

The CVD grown graphene used in this study have a significant electronic disorder which 
induces strong electron localization under IK (see characterization of bare graphene samples 
in Supplementary Figure SI and S3 and Supplementary Table SI). This results in a comple- 
tely different behavior : upon gating the tin decorated sample, one observes a superconductor 
to insulator transition (SIT). The SIT in 2D films of superconductors has been a very active 
field of research in the condensed matter community for the last twenty years [10J. It is of 
interest for the study of quantum phase transitions [TT] . but also to understand how su- 
perconductivity emerges in high-T c superconductors [12] . Recent advances in this field have 
involved the use of electric field to tune the SIT at constant disorder p^VTH] . 

The samples, presented in FigjTJ were fabricated by connecting CVD-grown graphene 
flakes transferred on oxidized silicon [15] (see Methods). A nominal thickness of 10 nm 
of tin was evaporated on the whole sample by thermal evaporation. Tin's dewetting (see 
FigjTJ) produces a self-assembled, non-percolating array of pancake-like nanoparticles [9] 
(see FigjTJ, b and d). The typical lateral size of a tin island is 80 nm, with a 13 nm gap 
in-between islands. 

Samples have shown two kinds of behavior depending on their room temperature sheet 
resistance. Devices having the lowest sheet resistance (typically < 10 kf2/D) showed super- 
conductivity for all gate voltages with a gate tunable transition temperature, similarly to 
what was reported for exfoliated samples [9]. Here, we will focus on the other type of devices, 
which exhibited a high sheet resistance (> 15 kfi/D) at room temperature. 

Upon cooling from room temperature down to 6K, the resistivity increase ranges from 
20% up to 100% at the charge neutrality point. This behavior is consistent with the enhanced 
weak localization and electron-electron interactions expected in 2D metals. Just above the 
critical temperature of bulk tin (T c 5ra = 3.7 K, black dotted line in Fig|2^), a 10% resistance 
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drop is observed (see Fig|2p), arising from superconducting fluctuations in graphene near the 
tin islands. Then a broad transition takes place, either towards a superconducting state at 
high electrostatic doping (for a voltage offset from the charge neutrality point AV > 45 V), 
or towards an insulating state for voltages closer to the charge neutrality point (AV < 10 V). 
In between these two gate voltages, the resistance levels off at low temperatures (Figj^), 
suggesting an intermediate metallic behavior. On the superconducting side, the system fol- 
lows an "inverse Arrhenius law" R oc R exp [Jf^)i as already reported in quench condensed 
granular films [16]. Significant fluctuations and leveling in the region just before supercon- 
ductivity sets in (see Figj3j bump near +10 V) are indicative of a percolative behavior. The 
amplitude of the critical current also supports this picture, as it was repetitively measured 
to be on the order of 1 //A, which corresponds to a critical current density of 5.10 -4 A/m, a 
value 2,000 times smaller than the value found in samples made from exfoliated graphene [9] 
or clean CVD graphene. 

Above 2K, R(T) curves at all gate voltages behave qualitatively the same. Reducing the 
temperature further, the curves near the charge neutrality point then reach a minimum 
at a gate-dependent temperature (red dotted line in Figj^), below which they start to 
increase sharply. This re-entrant behavior is reminiscent of what was observed in granular 
superconductors [T7] or in Josephson junctions arrays [18]. Our system is indeed similar 
to a granular superconductor in which the role of the intergranular media is taken up by 
graphene. In such systems, the SIT is driven by the competition between the charging 
energy E c of a superconducting island and the Josephson energy Ej = \-^-^~ tanh ' A 



2k B T J ' 

where Rn is the normal state resistance of the junction. Dissipative degrees of freedom, 
such as quasiparticles [19] or capacitive coupling to a 2DEG [20J like graphene [21J, lead to 
renormalization of the charging energy. Here, the leveling of resistance in the intermediate 
'metallic' regime is indicative of such dissipative processes [22J. Dissipation strength scales as 
R~x i an d Rn in turn depends on gate voltage and temperature, as can be seen by applying 
a magnetic field above the critical field (see Supplementary Figure S2). R N thus tunes 
both energies (Ej and E c ) simultaneously, unlike previously considered situations [20] |2T]. 
Despite this complex dependence, the phase diagram of the SIT (Fig. [2]i) shows that the 
boundary of the insulating region (red dotted line in Fig|2]D) can be related to a constant 
value of Rn- Interestingly, when measuring resistance at constant temperature as a function 
of gate(Fig.??), one sees a crossing point at V G = —20 V where the sheet resistance is 



of the order of j^p.-, the pair quantum of resistance. The universal value of the critical 
resistance at the transition was predicted by the so-called "dirty bosons model" [23] . Around 
this transition, the resistance varies by more than 7 orders of magnitude over a gate range of 
40 V (corresponding to a carrier density change of 3x 10 12 cm -2 ). This electrostatically driven 
transition shows a strongly insulating state, with exponential divergence of the resistance. 

The magnetoresistance curves in the insulating and superconducting regions are presented 
in Fig]4^L and show a peak, both in the superconducting and in the insulating regions. 
Such a non-monotonic behavior has been widely reported in superconducting thin films. By 
gating the sample, we observed a continuous crossover between different magnetoresistance 
regimes [24] . In the superconducting state (red curve in Figji]), the small resistance overshoot 
at intermediate magnetic field can be understood in terms of Galitski-Larkin correction to 
the conductivity [25]. The inflection point corresponds to the critical field expected in tin 
nanoparticles [26]. However, the behaviour in the insulating region (black curve in Fig|4^) 
cannot be explained with perturbation theory. Here, the resistance at intermediate field 
[B = 0.15 T) is about 40 times higher than the resistance in the normal state (B = 1 T). Such 
huge effects have been reported in amorphous thin films [271 EE] , and have been explained [29] 
to stem from the underlying nature of superconductivity in amorphous thin films, which is 
inhomogeneous [301 El] near the transition. In a granular system, since grains of different 
sizes have different critical fields, there exists an intermediate field where half of the grains are 
superconducting (S) and the other half are normal (N). N-S junctions can prevent percolation 
as they provide barriers for both quasiparticles and Cooper pairs. As we move away from 
the insulating region (by increasing Ay), the resistivity maximum is shifted towards higher 
magnetic fields. This indicates that islands with the smaller critical fields (Be ~ 0.15 T) 
play a crucial role in the percolation process, whereas deep in the superconducting region, 
coupling is established directly between other islands. 

Finally, the temperature dependance of the system when biased at AV = and for 
B = 0.15 T (i.e. in the region where Cooper pairs are localized) does not quite follow the 
activation law predicted by Beloborodov et al. |29| . Instead it shows an Efros-Schklovsky 
like behaviour T oc exp (p-) 2 with an activation energy T\ = 32.6 K (see Fig, 4^), suggesting 
that Coulomb interactions may play an important role in the transport. 

Going back to the gate-induced SIT, Figj5] shows how the data can be interpreted as a 
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quantum phase transition [TT], using finite-size scaling. The voltage range around the cri- 
tical gate voltage Vq c = -20 V was chosen the largest possible (±6V) while still retaining 
a universal exponent on both (insulating and superconducting) sides. Below 600 mK, the 
field effect at the transition shows a significant shift from the universal behavior. The pre- 
sence of puddles of normal electrons in the graphene sheet, which are sources of gap less 
excitations |21j, provides a dissipative environment. The system becomes more coupled at 
low temperatures to dissipative degrees of freedom, leading to a breakdown of the universal 
behavior as already observed in experiments involving MoGe films [22]. However, above 600 
mK, the critical resistance lies very close to the value predicted by the "dirty bosons" model : 
^ = 1.2, Rq = which is an indication that we are in the regime of low dissipation, 
where the dynamical critical exponent z is still equal to 1 [32]. The exponent zv has been 
evaluated using the two methods described in [33]. The first method is to multiply each 
curve by the factor t yielding the best collapse to the first curve, then fitting t to T~ zu . The 



second method is to take the slope of log 
in FigjBT) leads to zv = 1.05 ± 0.10 whi 



vs log (T 1 ). The first method (shown 



e the second method gives a value of zv = 1.18 
± 0.02 (see Supplementary Figure S4). The second method is probably more accurate, as 
it is based on data which are taken within the critical region. This value is close to other 
reported values of zv in thickness-tuned transition in Bi [33] . Note that this exponent differs 
from the expected exponent for classical percolation in 2D (zv=4/3). It is instead in good 
agreement with recent theoretical developments on the superfluid transition in disordered 
two-dimensional bosonic systems |34j . which can be understood as the percolation of phase 
coherent domains into a macroscopic superfluid. 



Unlike previously reported gate-induced SITs which either showed a partial SIT transition 
towards a weakly localized metal [121 H21 E2] or involved ionic gating, which freezes at low 
temperatures [H], tin-decorated CVD graphene can be gated continuously at low tempera- 
tures with a very strong transconductance (Fig (3]). This could have application, for example 
in transition-edge particle detectors. The recently demonstrated metal-insulator transition 
in ultra-clean graphene samples [4] also opens exciting new perspectives to probe the SIT in 
the opposite limit of very low disorder. The present experiment paves the way to the reali- 
zation of more complex graphene-based hybrid materials where graphene acts as a tunable 
medium or adjustable environment that controls the establishment of long-range electronic 
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orders, such as superconductivity or ferromagnetism. This experiment sheds light on the 
emergence of superconductivity in inhomogeneous superconductors, and demonstrates the 
potential of graphene as a versatile substrate for the realization of hybrid superconducting 
materials. 

METHODS 

Sample Preparation 

Graphene is grown using a chemical vapor deposition (CVD) technique on copper foils 
(typically 25 /xm thick, from Alfa-Aesar) following the methods described in [15] . During the 
growth, a flow of methane (CH4) provides the carbon feedstock, while forming gas (H^/Ar 
1:9) limits the reaction and only a single layer of graphene is obtained. After growth, the 
graphene is protected with a support layer of 1 /zm-thick polymethylmethacrylate (PMMA), 
and copper is etched away using a solution of 0.2 g/ml sodium persulfate (Na 2 S 2 8 ). The 
graphene remains attached to PMMA and floats in the solution. It is then carefully trans- 
ferred onto a wafer of degenerately doped oxidized silicon and PMMA is removed using 
acetone wash followed by thermal annealing at 380°C/lh under Argon atmosphere. Tin is 
deposited on the whole sample using room temperature Joule evaporation. Pd/Au electrodes 
are subsequently deposited using a millimeter-scaled metal foil stencil mask in a four-probe 
geometry aligned on top of the graphene sheet. Supplementary Figure SI shows a typical 
sample after fabrication. The fabricated samples were about 5 mm in length and 3 mm in 
width, sizes that could not previously be obtained using exfoliated graphene. Such a macro- 
scopic sample allowed us to get mesoscopic effects such as universal conductance fluctuations 
to be averaged out, which is crucial when studying how the phase transition scales. 

Several samples were measured, and the number of graphene layers was varied (from 1 
to 3), as well as the thickness of tin (8-20 nm). However, we did not see a direct correlation 
between these parameters and the behavior of the device. We could only relate it to the 
normal state sheet resistance of graphene. Only graphene showing a high sheet resistance 
at room temperature (>15 kf2) would behave as an insulator below the tin's supercon- 
ducting transition temperature (3.7 K). The other samples behaved like the ones studied by 
Kessler et al. [9], showing much higher critical current density and a gate-tunable Berezinsky- 
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Kozterlitz-Thouless transition towards a superconducting state at all gate voltages. 



Measurement setup 

The sample was thermally anchored to the mixing chamber of a 3 He/ 4 He dilution cryostat 
and connected to highly filtered measurement lines. The setup allowed the temperature to 
be continuously varied between 10 K and 0.03K. The differential resistance was recorded 
using a lock-in amplifier operated in a four-probe configuration at frequencies between 9 Hz 
and 37 Hz, with an excitation current of 1 nA. In the high impedance state, a two-probe, 
voltage-biased configuration was used using a Keithley 6430 electrometer or a Femto current- 
to- voltage converter to record the current. 
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Figure 1. Self-assembled graphene/tin nanohybrids. a : Atomic Force Micrograph of an 
l/j,m 2 area of the device, obtained by dewetting of an evaporated tin film of nominal thickness 
lOnm. b : Transmission Electron Micrograph (TEM) of a the decorated sample showing the self- 
organized network of tin nanoparticles (scale bar 200 nm) separated by clean graphene. c : TEM of 
the sample transferred on a membrane and observed at a grazing angle (the dashed line corresponds 
to the graphene surface) . The negative wetting angle of tin nanoparticles on graphene (black arrow) 
can be clearly seen, d Scanning Electron Micrograph (SEM) of the Sn nanoparticles network on 
graphene (scale bar is 300 nm). e,f Sketches of the device. Changing the gate voltage modulates 
the extension of phase coherent domains in graphene. g Photograph of the studied device. The 
dark region between the four electrodes is the decorated graphene sheet. The blue cast is due to 
the presence of tin nanoparticles on the whole surface. The enhanced contrast of the graphene 
sheet with respect to the silica sides comes from the difference of grain sizes and gaps between 
nanoparticles on graphene and on Si02- 
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Figure 2. Sheet resistance as a function of temperature for different gate voltages, a : 

Sheet resistance as a function of sample temperature for different gate voltages, plotted on a log 
scale. From top to bottom, voltage offsets from the charge neutrality point (Vr> = -36 V, see Figj3]) 
are AV = Vg - V D = 0, 3, 6, 8, 10, 12, 14, 16, 17, 18, 19, 20, 21, 22, 26, 31, 36, 41, 46, 56, 66, 76, 
86, 96 V. b : Same data (lower part) plotted on a linear scale to emphasize the behaviour between 
1 K and 4 K. The black dashed line indicates the critical temperature of tin. The red dashed line 
is a guide to the eye showing the minimum resistivity, c : Zoom on the critical region, d : Phase 
diagram of the superconducting to insulator transition, log (§^) vs gate voltage and temperature. 
I stands for "Insulating", and S for "Superconducting". We call superconducting the region where R 
is below the noise floor (R$ = 0.517). The black line is the iso- value of the normal state resistivity 
obtained at B = 1 T (see Supplementary Figure S2) which appeared the closest to the border 
between the two states. 
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Figure 3. Sheet resistance as a function of gate voltage for different temperatures. 

Gate voltage dependance of sheet resistance for the lowest temperatures. The vertical dotted line 
indicated the charge neutrality point. The horizontal one indicated the quantum of resistance for 
Cooper pairs i?n = (Top axis : carrier density calculated using the gate capacitance per unit 
area for 285 nm of Si0 2 C bg = 121 ^F/m 2 ). 



14 




Figure 4. Localization of Cooper pairs under a magnetic field, a : Sheet resistance as a 
function of magnetic field, measured at T = 300 mK. Black curve shows the magnetoresistance at 
the charge neutrality point. The red curve has been measured deep in the superconducting region, 
b : Temperature dependance of the sheet resistance at three different magnetic fields (indicated 
in the upper panel). The black lines are fits to the Efros-Schklovsky law, giving the following 
activation temperatures : T x = 7.8 K, T{ = 32.6 K, and Tf = 2.5 K for B = T, 0.15 T and 1 T, 
respectively. 
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Figure 5. Universal scaling of the transition Finite-size scaling analysis of the quantum phase 
transition, using data of Figj3j The best collapse was found using vz = 1.05 ± 0.05. Inset : subset 
of the data used for the scaling : 0.6 K < T < 1.3 K and a gate voltage interval of ± 6 V around 
the critical point Vq c = -20 V. 
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